function u = trapezoidal(u,N,x_max,k,num)
for i = 1:num
    u = trapezoidal_onestep(u,N,x_max,k);
end